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Abstract 

The critical thermodynamics of the two-dimensional N- vector cubic and MN 
models is studied within the field-theoretical renormalization-group (RG) ap- 
proach. The (3 functions and critical exponents are calculated in the five-loop 
approximation and the RG series obtained are resummed using the Borel- 
Leroy transformation combined with the generalized Pade approximant and 
conformal mapping techniques. For the cubic model, the RG flows for various 
iV are investigated. For N = 2 it is found that the continuous line of fixed 
points running from the XY fixed point to the Ising one is well reproduced 
by the resummed RG series and an account for the five-loop terms makes 
the lines of zeros of both (3 functions closer to each another. For the cubic 
model with N > 3, the five-loop contributions are shown to shift the cubic 
fixed point, given by the four-loop approximation, towards the Ising fixed 
point. This confirms the idea that the existence of the cubic fixed point in 
two dimensions under N > 2 is an artifact of the perturbative analysis. For 
the quenched dilute O(M) models (MN models with N = 0) the results are 
compatible with a stable pure fixed point for M > 1. For the MN model 
with M, N > 2 all the non-perturbative results are reproduced. In addition a 
new stable fixed point is found for moderate values of M and N. 

PACS Numbers: 75.10.Hk, 05.70.Jk, 64.60.Fr, ll.10.Kk 
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I. INTRODUCTION 



The two-dimensional model with iV-vector order parameter and cubic anisotropy is 
known to have a rich phase diagram; it contains, under different values of N and of the 
anisotropy parameter, the Ising-like and Kosterlitz-Thouless critical points, lines of the 
first-order phase transitions, and the line of the second-order transitions with continuously 
varying critical exponents (see, e. g. [1-3] for review). This model is related to many other 
familiar models in various particular cases, while for iV — > it describes the critical behavior 
of two-dimensional weakly disordered Ising systems. Moreover, exact solutions are known for 
the two-dimensional cubic model in the several limits such as an Ising decoupled limit, the 
limit of extremely strong anisotropy for N > 2 [2,4] and the replica limit N —> [2,5,6]. The 
correspondence, in particular regions of the phase diagram, with the iV-color Ashkin- Teller 
models, discrete cubic models, and planar model with fourth order anisotropy give further 
informations about the critical behavior. All these issues are reviewed in Ref. [3] and we 
will not repeat them here. These features make the two-dimensional iV-vector cubic model 
a convenient and, perhaps, unique testbed for evaluation of the analytical and numerical 
power of perturbative methods widely used nowadays in the theory of critical phenomena. 
The field-theoretical renormalizat ion-group (RG) approach in physical dimensions is among 
of them. 

Recently, the critical behavior of the two-dimensional iV-vector cubic model was explored 
using the renormalization-group technique in the space of fixed dimensionality [3]. The four- 
loop expansions for the /^-functions and critical exponents were calculated and analyzed using 
the Borel transformation combined with the conformal mapping and Pade-approximant 
techniques as a tool for resummation of the divergent RG series. The most part of predictions 
obtained within the renormalization group approach turned out to be in accord with known 
exact results. At the same time, some findings were quite new. In particular, for N > 2 
the resummed four-loop RG expansions for /3-functions were found to yield a cubic fixed 
point with (almost) marginal stability; this point does not correspond to any of the critical 
asymptotics revealed by exact methods ever applied. Although the stability properties of 
the cubic fixed point look very similar to those of its Ising counterpart, these points were 
found to lie too far from each other (for moderate N) to consider the distance between them 
as a splitting caused by the limited accuracy of the RG approximation employed. 

It is worthy to note that this situation is quite different from what we have in three 
dimensions. Indeed, for the three-dimensional cubic model the structure of the RG flow 
diagram is known today with a rather high accuracy. Recent five-loop [7-9] and six-loop 
[10,11] RG calculations certainly confirmed that for N > 2 the cubic fixed point does 
not merge with any other fixed point and governs the specific anisotropic mode of critical 
behavior, distinguishable from the Ising and Heisenberg modes (see, e. g. [12]). 

It is very desirable, therefore, to clear up to what extent the location of the cubic fixed 
point in two dimensions is sensitive to the order of the RG calculations and, more generally, 
whether this point really exists at the flow diagram or its appearance is the approxima- 
tion artifact caused by the finiteness of the perturbative series and by an ignorance of the 
confluent singularities significant in two dimensions [13-16]. 

Of prime interest is also the situation with the line of fixed points that should run, under 
N = 2, from the Ising fixed point to the XY one. Within the four-loop approximation, the 
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zeros of /^-functions for the 0(iV)-symmetric and anisotropic coupling constants form two 
lines that for N = 2 are practically parallel to each other and separated by the distance that 
is smaller than the error bar appropriate to the working approximation [3] . Will the higher- 
order contributions keep these two lines parallel? Will an account for the higher-order terms 
further diminish the distance between these lines or their splitting should be attributed, at 
least partially, to the influence of the singular terms just mentioned? 

To answer the above questions, it is necessary to analyze the critical behavior of the two- 
dimensional cubic model in the higher perturbative orders. Recently, the renormalization- 
group expansions for the two-dimensional 0(AQ -symmetric model were obtained within the 
five-loop approximation [15]. In the course of this study, all the integrals corresponding 
to the five-loop four-leg and two-leg Feynman graphs have been evaluated. This makes it 
possible to investigate the critical thermodynamics of anisotropic two-dimensional models 
with several couplings in the five-loop approximation. In this paper, such an investigation 
will be carried out for the two-dimensional iV-vector model with cubic anisotropy. 

For studying the effect of cubic anisotropies one usually considers the 4 theory [17,18]: 

r [ 1 N 1 N 1 

« = / d " x o £ [< W + r #] + 7i E («o + Mi) M , (1.1) 

[ i=l ' i,j=l ) 

in which the added cubic term breaks explicitly the O(N) invariance leaving a residual 
discrete cubic symmetry given by the reflections and permutations of the field components. 
This term favors the spin orientations towards the faces or the corners of an N- dimensional 
hypercube for v o < or v o > respectively. In two dimensions the effect of anisotropy is 
particularly important: systems possessing continuous symmetry do not exhibit conventional 
long-range order at finite temperature, while models with discrete symmetry do undergo 
phase transitions into conventionally ordered phase. 

In general, the model (1.1) has four fixed points: the trivial Gaussian one, the Ising 
one in which the N components of the field decouple, the ©-(AT) -symmetric and the cubic 
fixed points. The Gaussian fixed point is always unstable, and so is the Ising fixed point for 
d > 2 [17]. Indeed, in the latter case, it is natural to interpret Eq. (1.1) as the Hamiltonian 
of N Ising-like systems coupled by the 0(iV)-symmetric term. But this interaction is the 
sum of the products of the energy operators of the different Ising systems. Therefore, at the 
Ising fixed point, the crossover exponent associated with the 0(N) -symmetric quartic term 
is given by the specific-heat critical exponent aj of the Ising model, independently of N (for 
N = this argument is equivalent to the Harris criterion [19]). Since cc/ is positive for all 
d > 2 the Ising fixed point is unstable. Obviously, in two dimensions this argument only 
tells us that the crossover exponent at this fixed point vanishes. Higher order corrections 
to RG equations may lead either to a marginally stable fixed point [20] or to a line of fixed 
points. It was argued that for N > 3 the former possibility is realized, while for N = 2 the 
latter one holds (see Ref. [1-3] and references therein). 

The stability properties of the fixed points depend on N. For sufficiently small values of 
N, N < N c , the 0(iV)-symmetric fixed point is stable and the cubic one is unstable. For 
N > N c and d > 2, the opposite is true: the renormalization-group flow is driven towards 
the cubic fixed point, which now describes the generic critical behavior of the system. At 
N = N c , the two fixed points should coincide for d > 2. At d — 2, it is expected that 
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N c — 2 [21] and a line of fixed points connecting the Ising and the 0(2)-symmetric fixed 
points exists [1,3]. 

A generalization of the Hamiltonian Eq. (1.1) is obtained by considering N coupled 
O(M) vector models, instead of N Ising models. The resulting Hamiltonian (defining the 
so called MN model) is [17,18] 



where a,b = 1, . . . M and i, j = 1, . . . N. The continuous O(MN) symmetry is explicitly 
broken by the v term to Cjy x O(M) where is the discrete group of permutations 
of N elements. The presence of such discrete symmetry allows for a finite temperature 
phase transition under general values of N and M (differently from models with vq = 
0). For M — 1 it reduces to the Hamiltonian of the cubic model, but it has physical 
applications for M ^ 1 too. For iV — > and generic M, under the condition u < 
and Vq > 0, it describes the critical behavior of quenched dilute O(M) models [18,22]. 
For M = 2 and N = 2 (N = 3) and vq > 0, it is relevant for the second-order phase 
transition in planar (isotropic) antiferromagnets with complicated ordering as the three 
dimensional sinusoidal magnets TbAu2, DyC2, type-II antiferromagnets TbAs, TbP (type- 
Ill antiferromagnets K 2 IrCl 6 , sinusoidal TbD 2 ) and many others [23]. For M = N = 2 and 
v o < an exact mapping [24] relates the Hamiltonian (1.2) with the 0(2) x 0(2) symmetric 
one, describing the critical behavior of frustrated antiferromagnets with non-collinear order 
[25,18]. For N — > oo it describes O(M) models with constrains [26] (as O(M) models with 
annealed disorder). 

The Hamiltonian (1.2) has been largely studied in the framework of e expansion (e = 4— d) 
[27,28] and directly in three dimensions [28-32], mainly to understand the critical behavior of 
the antiferromagnets quoted above. Much less studies were devoted to the two dimensional 
case for M/ 1, since its main features can be understood by non-perturbative arguments. 
In fact, three fixed points always exist in the RG flow for any N and M: the always unstable 
Gaussian (u = v = 0), the O(M) (with u = 0) and the O(NM) (with v = 0). Their 
stability properties are known from exact arguments. As for the Ising fixed point in the 
cubic model, the stability of the O(M) fixed point is governed by the specific-heat critical 
exponent ckm of the O(M) model. Being (in two dimensions) au negative for all M > 1, 
the O(M) fixed point is stable for physical values of M. An exact argument ensures that 
the O(NM) fixed point is unstable when NM > N c = 2 [21,33], thus it is always unstable 
for physical values (apart from the replica limit iV — > 0, but in this case it is not reachable 
from the physical initial condition with Uq > 0). In the framework of e expansion another 
fixed point (called mixed) exists [17]. It is the generalization to M ^ 1 of the cubic one. 
Its analytic continuation to d — 2 is expected to be in the region u < and v o > for 
N, M > 2 and in the uo, fo > region (unphysical) in the replica limit. It is expected to be 
unstable in both the cases. 

Anyway, these arguments still do not completely characterize the RG flow of the MN- 
model. Indeed for N, M > 2 and Vq < the situation is more controversial. The e expansion 
indicates that no fixed point exists and consequently it was longly believed that no continuous 
transition can take place in these systems. Nevertheless, for iV = M = 2 the quoted mapping 
between the MN and the 0(2) x 0(2) models shows the presence of a fixed point (found in 
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(1.2) 
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the 0(2) x 0(2) model at four [34] and five loops [35,36]) in the region vq < 0. Consequently 
the transition should be second order in the chiral universality class. The natural question 
arising is whether the fixed point found at N = M = 2 is a peculiar feature or it persists to 
some larger values of N and M (not too large, since for N >> 1, at fixed M, we expect that 
the e expansion can be smoothly continued up to d — 2). An answer to this question can be 
found within higher order perturbative expansion. 

The paper is organized as follows. In Section II the five-loop contributions to the 
renormalizat ion-group functions are calculated and the singularities of Borel transforms 
of the renormalization-group series are discussed for the general MN model. Section III is 
devoted to the analysis of the critical behavior of the model in the five-loop approximation. 
The existence of a cubic fixed point for N > 3 is investigated. The critical behavior of 
quenched random O(M) models is considered. Attention is paid to the reproducibility of a 
continuous line of fixed points, predicted earlier for the planar (N — 2, M — 1) model, in 
the framework of the field-theoretical RG approach. The critical exponents along this line 
are evaluated in this section as well. Finally, the critical behavior for N, M > 2 is addressed. 
In Section IV we summarize the main results obtained and make concluding remarks. 



II. FIXED DIMENSION PERTURBATIVE EXPANSIONS 



A. Renormalization of the theory 



The fixed-dimension field-theoretical approach represents a powerful procedure in the 
study of the critical properties of three-dimensional systems belonging to the O(N) and 
more complicated universality classes (see, e.g., Ref. [18,37,38]). In this approach one per- 
forms an expansion in powers of appropriately defined zero-momentum quartic couplings 
and renormalizes the theory by a set of zero-momentum conditions for the (one-particle 
irreducible) two-point and four-point correlation functions. For the MN model they read: 



rg«(p) = ^I 1 m 2 +p 2 + 0(p 4 ) 



(2.1) 



where 8 aijbj = 5 ab 5ij, and 



2^2 



U 



— {&ai,bjbck,dl + $ai,ckhj,dl + 8 a i,dl$bj,ck) + 



- 5ij5 ik 5ii(5 ab 5 cd + 5 ac 5 bd + 5 ad 5 bc ) 



(2.2) 



They relate the inverse correlation length (mass) m and the zero-momentum quartic cou- 
plings u and v to the corresponding Hamiltonian parameters r, tt , and Vq: 

u q = m 2 uZ u Z^ 2 , v = m 2 vZ v Z^ 2 . (2.3) 

In addition, one introduces the function Z t defined by the relation 

T a ]f j (0) = S ai , bj Zi 1 , (2.4) 

where r^ 1,2 ) is the (one-particle irreducible) two-point function with an insertion of |0 2 . 
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From the perturbative expansion of the correlation functions T^ 2 \ , and I^ 1,2 ) and the 
above relations, one derives the functions Z^Ujv), Z u (u,v), Z v (u,v), and Z t (u,v) as double 
expansions in u and v. 

The fixed points of the theory are given by the common zeros of the /3-functions 



P u {u,v) = m — 



UQ,V0 



P v {u,v) = m — 
dm 



(2.5) 



The stability properties of the fixed points are controlled by the eigenvalues u>i of the matrix 

fd(3 u (u,v) d(3 u (u,v)\ 

(2.6) 



du dv 
dfl v (u,v) d(3 v (u,v 



\ du 



dv ' 



computed at the given fixed point: a fixed point is stable if both eigenvalues are positive. 
The eigenvalues Ui are related to the leading scaling corrections, which vanish as ~ \t\ Ai 
where A, = vloi. If u.i has a non- vanishing imaginary part, the approaching to the FP is 
spiral-like and the FP is called focus [35]. 
One also introduces the functions 



dhiZ, 



d\nm 



d In Z^t d In Z^ t 



u ,v 



du 



dv 



in terms of that, the critical exponents are obtained as 

V = V<p(u*,v*), 

u=[2-r ]4> (u*,v*)+Vt(u*,v*)r 1 
7 = 1/(2-77), 



(2.7) 



(2.8) 
(2.9) 
(2.10) 



where (u*,v*) is the location of the stable fixed point. 



B. The five loop series 

We calculate the two-dimensional perturbative RG functions of the MN model in the 
zero-momentum massive renormalization approach introduced before up to five loops. Cor- 
responding contributions to the functions Eqs. (2.1), (2.2), and (2.4) are given by 162 
diagrams for the four-point correlators and by 26 graphs for the two-point one [39]. We 
handle them with a symbolic manipulation program, which generates the diagrams and 
computes the symmetry and group factors of each of them. The RG functions are written 
in terms of the rescaled couplings 

U = ~Y R nm % V = ^ R M v, (2.11) 

where R K = 9/(8 + K). 
The obtained series are 
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_ 2 2(M + 2) 



flu = -U + U + 



(M + 



-uv 



i+j>2 



(2.12) 



= -u + v 2 + 



T7V7 UV + V ^ bi^UV 3 , 

8 + NM ^ 3 



(2.13) 



i+j>2 



eWvh? 



[2.U) 



where 



2(2 + iVM)_ 2 + M_ „ „ 
m ~ ~ (8 + NM) U ~ 8TM" + j+ % 2 6ij 



(2.15) 



16n 



16tt 



(2.16) 



The coefficients b\j\ and efj are reported in the Tables I, II, III, and IV. Note 
that due to the rescaling (2.16) and (2.11), the matrix element of fl are the double of the 
derivative of /3 with respect to u and v. 

We have verified the exactness of our series by the following relations: 

(i) j3u(% 0), rf^(u, 0) and r) t (%0) reproduce the corresponding functions of the O(NM)- 
symmetric model [15,40]. 

(ii) j3v(0,v), 7)^(0, v) and r) t (0, v) reproduce the corresponding functions of the O(M)- 
symmetric 4 theory [15,40]. 

(iii) The following relations hold close to Heisenberg 0(NM) [41,42] 



du 



MN + 2 M + 8 driM 



du 



(«,0) 



(u,0) 



M + 2 MN + 8 dv 

(«,o) 



dv 



(u,0) 

MAT + 2 M + 8 dfc 



M + 2 MN + 8 dv 



(u,0) 



and to the O(M) fixed points [42] 



dv 



dh 



(0,v) 



du 



MN + 8 dfc 



dv 



(0,v) 



(0,v) 



M + 8 du 
MN + 8 dricp 



M + 8 du 



(0,v) 



(0,5) 



(iv) The following relations hold for N — 1 and generic M: 

Pu(u, x - u) + (3v(u, x - u) = j3y(0, x) , 
r),p(u,x - u) = 77^(0, a;), 
r] t (u,x-u) = r} t (0,x). 



(2.17) 
(2.18) 



(2.19) 
(2.20) 



(2.21) 
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(v) For N = 2 and M — 1, one easily obtains the identities [43,3] 



+ p, -v) + -fc(u + p, -v) = Pu(u, v) , (2.22) 

j3y(u + |u, -U) = -f3v(u,v), 
rj^u+p, -v) = rfo(u,v), 
T) t (u+ p, -v) = 7) t (u,v). 

These relations are also exactly satisfied by our five- loop series. Note that, since the Ising 
fixed point is (0, g}), and g\ is known with high precision [44] 

gj = 1.7543637(25), (2.23) 

the above symmetry gives us the location of the cubic fixed point: (\g}, —<?/)■ 

(vi) In the large-iV limit the critical exponents of the mixed fixed point are related to 
those of the 0(M) model: r\ = i] M and v = u M [26,45] (note that in general the latter 
equivalence holds only for M > 1 and d — 2, in the general case Fisher renormalization of 
exponents [46] has to be taken into account [45]). One can easily see that, for iV — > oo, 
r)<j,(u,v) = r)<j,(0, v), where 7^(0, t>) is the perturbative series that determines the exponent 
r/M of the 0(M) model. Therefore, the first relation is trivially true. On the other hand, 
the second relation v = um is not identically satisfied by the series, and is verified only at 
the critical point [10]. 

(vii) For M — 0, the iV-independent series satisfy [47] 

/3u(u, x - u) + j3v(u, x - u) = flu(x, 0) , 

7)4,(11, x-u) = 7),p(x,0) , 

r k (u,x-u)=r k (x,0). (2.24) 

(viii) For M — > oo the series, as expected, reduce to 

f3u = — u + u 2 + 2uv, j3v = —v + U 2 , 

7) t = -2u -v, 7)4 = 0, (2.25) 

with two Gaussian and two spherical (i.e. O(oo)) fixed points (as the O(N) x 0(M) model 
[48]). 

(ix) For N = M = 2 they agree, according to the exact mapping [24]: 

U = Uch + ~7~ j V = ~^Vch , (2.26) 

2 6 

with the five-loop series of the 0(2) x 0(2) model [36], written in terms of w c h and U c h- 

(x) The series reproduce the previous four- loop results for M — 1 [3]. 

C. Resummations of the series and analysis method 

The obtained RG series are asymptotic and some resummation procedure is needed 
in order to extract accurate numerical values for the physical quantities. Exploiting the 
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TABLE I. The coefficients b^' , cf. Eq. (2.12). 



2.0 -0.588581 - 0.127593 MN 

1.1 -0.621503 - 0.310751 M 
0,2 -0.144054 - 0.0720268 M 

3.0 0.719309 + 0.204598 MN + 0.00685909 M 2 N 2 

2.1 1.19188 + 0.595938M + 0.0598976MV + 0.0299488M 2 7V 

1.2 0.661318 + 0.410885M + 0.0401129M 2 + 0.00432058MN + 0.00216029M 2 AT 
0,3 0.131463 + 0.0886454 M + 0.0114569 M 2 

4.0 -1.157005 - 0.397981 MN - 0.0273887 M 2 N 2 - 0.0000135411 M 3 N 3 

3.1 -2.6265 - 1.31325 M - 0.26635 MN - 0.133175 M 2 N + 0.00011597 M 2 N 2 + 0.0000579851 M 3 N 2 

2.2 -2.31478 - 1.56359 M - 0.203103 M 2 - 0.0493674 MN - 0.0244594 M 2 N + 0.000112162 M 3 N 

1.3 -0.970048 - 0.704633 M - 0.10947 M 2 + 0.00016724 M 3 - 0.00381205 MN - 0.00238253 M 2 N - 0.000238253 M 3 N 
0,4 -0.166668 - 0.124257 M — 0.0206675 M 2 - 0.000103042 M 3 

5.0 2.26883 + 0.901262 MN + 0.0889695 M 2 N 2 + 0.00136005 M 3 N 3 - 6.9088710" s M 4 Af 4 

4.1 6.58636 + 3.29318 M + 0.970818 M N + 0.485409 M 2 N + 0.0167572 M 2 N 2 + 0.00837862 M 3 N 2 
+0.000014758 M 3 N 3 + 7.37899 10" 6 M 4 N 3 

3.2 8.02544 + 5.66644 M + 0.826859 M 2 + 0.359772 M N + 0.221805 M 2 N + 0.0209597 M 3 N + 0.00113766 M 2 N 2 
+0.000617202 M 3 N 2 + 0.0000241853 M 4 N 2 

2.3 5.23732 + 4.0716M + 0.770463M 2 + 0.0219952M 3 + 0.058081MTV + 0.039797M 2 A r + 0.0054257M 3 V + 0.000023694M 4 7V 

1.4 1.84567 + 1.47827 M + 0.302198 M 2 + 0.0122745 M 3 + 0.000017524 M 4 + 0.00491992 M N + 0.00347 M 2 N 
+0.000475809 M 3 N - 0.0000146041 M 4 N 

0,5 0.278724 + 0.226356 M + 0.0478936 M 2 + 0.00219097 M 3 - 3.68575 10 -6 M 4 



property of Borel summability of <ft 4 theories in two and three dimensions [49], we resum the 
divergent asymptotic series by a Borel transformation combined with an analytic extension 
of the Borel transform to the real positive axis. This extension can be obtained by a Pade 
approximant or by a conformal mapping [50] which maps the domain of analyticity of the 
Borel transform onto a circle (see Refs. [50,37] for details). 

The conformal mapping method takes advantage of the knowledge of the large order 
behavior of the perturbative series F(u,z) = J2k fk(z)u k [50,37] 

f k {z) ~ k\ (-a(z)) k k b [l + Oik' 1 )] with a(z) = -l/u h {z), (2.27) 

where Ub(z) is the singularity of the Borel transform closest to the origin at fixed z = v/u, 
given by [10,3,30] 

1 2N R 

— — =- a (R MN + R M z) for < z and z <- ~Er^, ( 2 -28) 

u b (z) N + 1 R M 

^ ( 7-> 1 n \ r r> RmN 

—a ( Rmn + jjRmz J tor > z > 



u h {z) ^ m " N m > N+l R 



M 



where a = 0.238659217 . . . [37]. Note that the condition of Borel summability (that coincides 
with the mean-field boundness condition) is z > z\ = —Rmn/Rm for v < 0. So in this region, 
even if the second of Eq. (2.28) takes into account the singularity of the Borel transform 
closest to the origin, there is another singularity on the real positive axis that makes the 
series not Borel summable. The situation is the same as for O(M) x O(N) under > 
[51,34,35]. As in the latter case we resum the series even where they are not Borel summable. 
Although in this case the sequence of approximation is only asymptotic, it should provide 
reasonable estimates as long as we are taking into account the leading large-order behavior 
(i.e. as long as z > z 2 = - n+i R rJ )- 
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TABLE II. The coefficients b^f , cf. Eq. (2.13). 



2.0 —1.14424 — 0.0720268 M AT 

1.1 —1.62169 — 0.310751 M 
0,2 -0.588581 - 0.127593 M 

3.0 1.68536 + 0.162556 M N — 0.00251242 M 2 N 2 

2.1 3.65454 + 0.776108 M + 0.036274 M N — 0.0011187 M 2 N 

1.2 2.74577 + 0.729746 M + 0.0159795 M 2 
0,3 0.719309 + 0.204598 M + 0.00685909 M 2 

4.0 -3.15852 — 0.410304 M N — 0.00388438 M 2 N 2 — 0.00012569 M 3 N 3 

3.1 -9.25529 - 2.09353 M - 0.273324 MAT- 0.0473612 M 2 N + 0.00103599 M 2 N 2 - 0.000247572M 3 A r2 

2.2 -10.5999 - 3.25027 M - 0.159554 M 2 - 0.0210851 M N — 0.00248042 M 2 N - 0.0000643977 M 3 N 

1.3 -5.61836 - 1.86565 M - 0.116276 M 2 + 0.0000353371 M 3 
0,4 -1.157 - 0.397981 M - 0.0273887 M 2 - 0.0000135411 M 3 

5.0 7.02662 + 1.14357 M N + 0.0316503 M 2 N 2 - 0.000233057 M 3 N 3 - 7.79352 10~ 6 M 4 N 4 

4.1 26.007 + 6.15895 M + 1.33378 M N + 0.286892 M 2 N - 0.0021517 M 2 N 2 - 0.00120819 M 3 N 2 
+0.0000320327 M 3 N 3 - 0.0000252216 M 4 N 3 

3.2 40.1395 + 13.492 M + 0.87123 M 2 + 0.403601 M N + 0.0984158 M 2 N - 0.00021126 M 3 N — 0.000514422 M 2 N 2 
-8.56408 10~ 6 M 3 N 2 - 0.0000250754 M 4 N 2 

2.3 32.1867 + 12.0427M + 1.0329M 2 + 0.00813502 M 3 - O.OOegggeSMAT - 0.004025M 2 Af 
-0.000484565M 3 iV - 3.956 lO- 6 !^ 4 ^ 

1.4 13.3343 + 5.18122 M + 0.485923 M 2 + 0.00596934 M 3 + 3.27121 10" 6 M 4 
0,5 2.26883 + 0.901262 M + 0.0889695 M 2 + 0.00136005 M 3 - 6.90887 10 -8 M 4 



These results do not apply to the case N = 0. Indeed, in this case, additional singularities 
in the Borel transform are expected [52]. 

For each perturbative series R(u,v), we consider the following approximants [50] 

E(R) p (a,b;u,v) = £ B k (a,b;v/u) f°° dtt^ ^^J^" (2.29) 
t^o Jo [1 -y{ut;v/u)] a 

where 



\ 1 — x/u b (z) - 1 

y(*; *) = -i == ■ (2-30) 

yjl - X/U h {z) + 1 

The coefficients B k are determined by the condition that the expansion of E(R) p (a,b;u,v) 
in powers of u and v gives R(u, v) to order p. Within this method we end up with several 
values for each resummed quantity, depending upon the free parameters a and b. To have a 
single estimate we search for the values of a opt and 6 opt minimizing the difference between the 
highest orders and, as usual [50,10,3], we consider as final estimate (error bar) the average 
(the variance) of the approximants with a G [a opt — Act, a opt + Act] and b e [b opt — Ab, b opt + 
Ab], with Aa = 2 and Ab = 3 (see e.g. Ref. [10] for a discussion about the effectiveness of 
such choice of Aa and Ab). 

Within the second resummation procedure, the Borel-Leroy transform is analytically 
extended by means of a generalized Pade approximant, using the resolvent series trick (see, 
e. g. [9]). Explicitly, once introduced the resolvent series of the perturbative one R(u,v) 

n p 

P(R) (% v, b,\) = y\ n y . k ' n -u n - k v k , (2.31) 

1 ^ ^ r fr o r(n + 6 + l) 1 ; 
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TABLE III. The coefficients e\y , cf. Eq.(2.14). 




2,0 

1,1 
0,2 
3,0 
2,1 
1,2 
0,3 
4,0 
3,1 
2,2 
1,3 
0,4 
5,0 
4,1 



1,4 
0,5 



3,2 



2,3 



0.0226441 + 0.011322 M N 
0.0452882 + 0.0226441 M 
0.0226441 + 0.011322 M 

-0.00119855 - 0.000749094 M N — 0.0000749094 M 2 N 2 
-0.00359565 - 0.00179783 M - 0.000449457 M N — 0.000224728 M 2 N 
-0.00359565 - 0.00224728 M - 0.000224728 M 2 
-0.00119855 - 0.000749094 M - 0.0000749094 M 2 

0.00633062 + 0.00445834 M N + 0.000618261 M 2 N 2 - 0.0000141266 M 3 TV 3 

0.0253225 + 0.0126612 M + 0.00517211 M N + 0.00258606 M 2 N - 0.000113013 M 2 N 2 - 0.0000565063 M 3 N 2 
0.0379837 + 0.0254875 M + 0.00324782 M 2 + 0.00126253 M N + 0.000461744 M 2 N - 0.0000847595 M 3 N 
0.0253225 + 0.0178334 M + 0.00247304 M 2 - 0.0000565063 M 3 
0.00633062 + 0.00445834 M + 0.000618261 M 2 - 0.0000141266 M 3 

-0.00722953 - 0.00550947 M N — 0.000978196 M 2 N 2 - 0.0000178225 M 3 N 3 - 1.20103 10~ 6 M 4 N 4 
-0.0361476 - 0.0180738 M - 0.00947355 M N — 0.00473677 M 2 N - 0.000154204 M 2 N 2 - 0.0000771022 M 3 N 2 
-0.0000120103 M 3 N 3 - 6.00517 10~ 6 M 4 N 3 

-0.0722953 - 0.0502342 M - 0.0070433 M 2 - 0.00486051 M N — 0.00273867 M 2 N - 0.000154207 M 3 N+ 
5.95387 10^ 9 M 2 N 2 - 0.0000240177 M 3 N 2 - 0.0000120103 M 4 N 2 

-0.0722953 - 0.0544783 M - 0.00942071 M 2 - 0.000127693 M 3 - 0.000616447 M N — 0.000361246 M 2 N 
-0.0000505316 M 3 N - 0.0000120103 M 4 N 

-0.0361476 - 0.0275474 M - 0.00489098 M 2 - 0.0000891125 M 3 - 6.00517 10" 6 M 4 
-0.00722953 - 0.00550947 M - 0.000978196 M 2 - 0.0000178225 M 3 - 1.20103 10~ 6 M 4 



which is a series in powers of A with coefficients being uniform polynomials in u,v. The 
analytical continuation of the Borel transform is the Pade approximant [N/M] in A at A = I. 
Obviously, the sum for each perturbative series depends on the chosen Pade approximant 
and on the parameter b. We consider for the final estimates all the non-defective (i.e. having 
no singularities on the real positive axis) Pade approximants at the value of b minimizing 
the difference between them. 

An important issue in the fixed dimension approach to critical phenomena (and in general 
of all the field theoretical methods) concerns the analytic properties of the RG functions. 
As shown in Ref. [16] for the O(N) model, the presence of confluent singularities in the 
zero of the perturbative f3 function causes a slow convergence of the results given by the 
resummation of the perturbative series to the correct fixed point value. The O(N) two- 
dimensional field-theory estimates of physical quantities [50,15] are less accurate than the 
three-dimensional ones [50] due to the stronger non-analyticities at the fixed point [16,13,14], 
to say nothing about the stronger growth of the series coefficients themselves. In Ref. [16] 
it is shown that the non-analytic terms may cause large imprecision in the estimate of the 
exponent related to the leading correction to the scaling u. At the same time, the result for 
the fixed point location turns out to be a rather good approximation compared with those 
coming from different techniques. Non-analyticities cause slow convergence to the exact 
critical exponents too. 



We start the analysis of five-loop perturbative series in the case when one coordinate 
of the fixed point (FP) vanishes. Since in this case a lot of exact results are known, these 



III. ANALYSIS OF FIVE LOOP SERIES 



A. Stability of the O(NM) and 0{M) fixed points 
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TABLE IV. The coefficients e^' , cf. Eq.(2.15). 



2.0 0.166698 + 0.0833489 M N 

1.1 0.333395 + 0.166698 M 
0,2 0. 166698 + 0.0833489 M 

3.0 -0.132662 - 0.0893945 MN — 0.0115318 M 2 N 2 

2.1 -0.397986 - 0.198993 M - 0.0691908 M N - 0.0345954 M 2 N 

1.2 -0.397986 - 0.259542 M — 0.0302748 M 2 - 0.00864116 M N — 0.00432058 M 2 N 
0,3 -0.132662 - 0.0893945 M - 0.0115318 M 2 

4.0 0.172999 + 0.128715 MN + 0.0212857 M 2 N 2 + 0.0000889159 M 3 N 3 

3.1 0.691995 + 0.345998 M + 0.168863 M N + 0.0844315 M 2 N + 0.000711327 M 2 N 2 + 0.000355664 M 3 N 2 

2.2 1.03799 + 0.711189 M + 0.0960966 M 2 + 0.0611014 M N + 0.0316177 M 2 N + 0.000533495 M 3 N 

1.3 0.691995 + 0.507236 M + 0.0803778 M 2 - 0.000120843 M 3 + 0.0076241 M N + 0.00476506 M 2 N + 0.000476506 M 3 N 
0,4 0.172999 + 0.128715 M + 0.0212857 M 2 + 0.0000889159 M 3 

5.0 -0.285954 - 0.231865 MN - 0.0488718 M 2 N 2 - 0.0022088 M 3 N 3 + 2.48471 10~ 6 M 4 N 4 

4.1 -1.42977 - 0.714885 M - 0.444443 MN - 0.222221 M 2 N - 0.0221377 M 2 N 2 - 0.0110688 M 3 N 2 
+0.0000248471 M 3 V 3 + 0.0000124236 M 4 JV 3 

3.2 -2.85954 - 2.02478 M - 0.297507 M 2 - 0.293871 MN- 0.189024 M 2 N - 0.0210442 M 3 N 
-0.00218691 M 2 N 2 - 0.00104376 M 3 N 2 + 0.0000248471 M 4 N 2 

2.3 -2.85954 - 2.23824 M - 0.433491 M 2 - 0.0146269 M 3 - 0.0804108 MN - 0.055227 M 2 N 
-0.0074611 M 3 N + 0.0000248471 M 4 N 

1.4 -1.42977 - 1.14949 M — 0.237419 M 2 - 0.0100924 M 3 - 0.0000167847 M 4 - 0.00983986 M N - 0.00694 M 2 N 
-0.00095162 M 3 N + 0.0000292083 M 4 N 

0,5 -0.285954 - 0.231865 M - 0.0488718 M 2 - 0.0022088 M 3 + 2.48471 lO" 6 M 4 



calculations will be a check of the goodness of our resummation procedure and of the error 
made considering only five loops in the perturbative expansion. 

The location of the FP on the axis v = (that is equivalent to u — 0, with replacing 
NM — > M) was already studied in Ref. [15]. Also the critical exponents were considered 
there. The agreement of these results with the available exact ones was satisfactory (see also 
footnote [40]). What still remains to be computed is the stability of the FP's with respect 
to perpendicular perturbation. 

First of all, we analyze the stability properties of the 0(A r M)-symmetric fixed point. 
Since duPv\(u,o) = 0, the stability with respect to an anisotropic cubic perturbation is given 
by u v (u) = 2(9t7/3f7|( Si0 ), whose five-loop expansion is (n = NM) 

u v (u) 12 92.6834 + 5.83417n_ 2 1228.63 + 118.503n - 1.83156ra 2 _, 

= — 1 -f- 7/ — II A- II 

2 n + 8 (n + 8) 2 (n + 8) 3 

20723.1 + 2692.00n + 25.4854n 2 - 0.824655^3 

u 



(n + 8) 4 

U 5 . (3.1) 



414915.7 + 67526.8 n + 1868.92 n 2 - 13.7618 n 3 - 0.4602 n 4 _. 



(n + 8) 5 

This series must be calculated at the 0{n) FP located at u = u*. In Table V we report 
the results for uo v for several values of n, obtained resumming the series (3.1) at the FP's 
calculated in Ref. [15]. The 0(n) FP results unstable for n > 3. For n = 2 our result is 
compatible with the expected result uj v = 0, which is essential in the context of a continuous 
line of the fixed points. 

We can also use these results to discuss the nature of the multicritical point in two- 
dimensional models with symmetry 0(N 1 ) © O(A^) [53]. Indeed, they allow to exclude that 
the multicritical point has enlarged symmetry 0(N 1 + N 2 ) if Ni + N 2 = n > 2 [33]. In RG 
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TABLE V. Half of the exponent u> v at the 0(n) fixed point. CM is the value obtained using 
conformal mapping technique and PB the one using Pade-Borel. 



n 


CM 4-loop 


PB 4-loop 


CM5-loop 


PB 5-loop 


2 


0.03(3) 


0.06(4) 


0.025(40) 


0.00(5) 


3 


-0.08(3) 


-0.07(3) 


-0.10(6) 


-0.10(5) 


4 


-0.18(4) 


-0.17(5) 


-0.17(4) 


-0.20(4) 


8 


-0.45(5) 


-0.44(6) 


-0.48(4) 


-0.50(5) 



terms, this can generally occur if the 0(Aq + N 2 ) fixed point has only two relevant O(Aq)© 
0(A^2)-symmetric perturbations. But, when n > 2, the instability of the O(n) fixed point 
with respect to the cubic perturbation shows that at least one extra relevant perturbation 
exists (see for details Ref. [33]). 

Then we focus our attention on the stability properties of the 0(M) fixed point. Also in 
this case the stability is given by u u (v) = 2du/3u\(o,v), evaluated at the O(M) fixed point v*. 
As expected, the series uj u (v) is independent of N 

ojJv) „M + 2 5.83417(M + 2) , (47.9183 + 8.35207M)(M + 2) „ 

v ' - _ -L 9 77 - -77 4- - — -77 

2 M + 8 (M + 8) 2 (M + 8) 3 

(546.755 + 134.247M + 0.67 606M 2 )(M + 2) _ 4 



(M + 8) 4 

(8229.19 + 2568.45M + 129.81M 2 - 0.2176M 3 )(M + 2)_ 5 
+ (M + 8)* V ■ (3 ' 2) 

For M — 1, the fixed point value of this exponent is oo T u /2 = —0.09(8), using the conformal 
mapping method, and —0.08(10) with the Pade-Borel. We note that the [4/1] approximant 
with b = 1 leads to = —0.031. These values are compatible with the exact known 

result otj = 0. For M > 2, we always find uo u > 0, in agreement with the exact results 
predicting the stability of the 0(M) fixed point. 

An alternative approach to the analysis of the series Eqs (3.1) and (3.2) is the so called 
pseudo-e expansion [54]. Within this method, one multiplies the linear term in the j3 func- 
tions by a fictitious small parameter r, finds the common zeros of the f3 functions as series in 
r, and finally the exponents (as series in r) are obtained introducing such expansions in the 
corresponding RG functions. This approach has twofold advantage. First, in the estimates 
of the critical exponents the cumulation of errors due to not exact knowledge of the FP and 
that of the proper uncertainty of the RG function is avoided [50]. Second, it allows for very 
precise estimates of the marginal number of components since these are series in r that must 
be evaluated at r = 1 [32,31,55,56]. We use this method only for the second purpose. 

Imposing the vanishing of u v at the 0(NM) fixed point, the series for N c determining 
the relevance of cubic anisotropy (see the Introduction) is: 

N c = 4 - 2.2504r + 0.6230r 2 - 0.7509r 3 + 1.1761r 4 . (3.3) 

Analogously, imposing u u = at the O(M) FP leads to 

M c = 4 - 4.5008r + 2.1693r 2 - 1.2165r 3 + 1.3055r 4 , (3.4) 
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determining the relevance of an energy-like interaction, as e.g. the relevance of quenched 
randomness. Such series do not behave asymptotically with a factorial growth of the coef- 
ficients (at least up to the considered order), thus simple Pade approximants without Borel 
resummation are expected to give reliable estimates. This in fact turned out to be the case 
for their three dimensional analogs [32,31] and also for more complicated series as those for 
the marginal number of components of O(N) x O(M) [55] and U(N) x U(M) [56] models. 
The Pade table for the series Eq. (3.3) at r = 1 is 



4 
2.56 
2.32 
2.08 
2.21[3.i] 



1.75 
2.24 

5.76[o.9] 
2.17 



2.37 1.62 2.53 
2.03 2.08 
2.08 



(3.5) 



where we indicated in small brackets the closest pole to the origin on the real axis (when 
exists). Whenever this pole is close to r = 1, the approximant has not to be considered in 
the average procedure. This Pade table leads to the estimates N c = 2.10(7) that includes 
all the four- and five-loop estimates without poles (excluding the first line and column that 
are unreliable). This is in quite good agreement with the known exact result N c = 2 [21], 
signaling about the high effectiveness of the perturbative expansion technique at the five-loop 
level. 

Similarly, the Pade table for the series Eq. (3.4) is 



4 -0.5 1.66 0.45 1.75 
1.88 0.96 0.89 1.08 

1.40 0.88 [8 .2] 0.95 
1.19 1.09 
1.13 



(3.6) 



The same averaging procedure as before leads to M c = 0.99(10), that perfectly agrees with 
the exact value M c = 1 at which the specific-heat exponent is vanishing. This demonstrates 
under u = good approximating properties of the five-loop series. 

We analyzed the series (3.3) and (3.4) with the Pade-Borel resummation as well, having 
obtained equivalent results. So, indeed at the considered order these series do not behave 
as asymptotic. 



B. The cubic model for N > 3 

In this Section we consider the existence of the cubic fixed point for N > 3, previously 
found in the four-loop approximation [3]. It was claimed in Ref. [3] that the quite peculiar 
features of this fixed point (like the marginal instability) make its existence quite doubtful 
and that it might be an artifact of the relatively short series available at that time. Now we 
are in a position to analyze longer series and to further confirm or to reject this statement. 

The results obtained with the conformal mapping methods are reported in Table VI 
together with the previous four loop results [3], in order to make the comparison visible at 
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TABLE VI. Location of the apparent cubic fixed point for some N > 3. 



N 


CM 4-loop 


CM 5-loop 


Pade [4/1] b = 1 


3 


[0.83(12), 1.12(9)] 


[0.54(6), 1.35(4)] 


[0.050,1.757] 


4 


[0.54(10), 1.43(8)] 


[0.32(5), 1.58(4)] 


[0.031,1.774] 


8 


[0.24(8), 1.72(10)] 


[0.14(4), 1.74(4)] 


[0.015,1.788] 



first sight. This Table shows that the cubic fixed point drastically moves towards the Ising 
fixed point with increasing the order of perturbation theory from four to five loops. It may 
be considered as an argument in favor of the statement that, within the exact theory, the 
cubic and the Ising fixed points coincide. On the other hand, both at four and at five loops, 
the quoted errors are less than the difference between the two estimates. This leads to the 
conclusion that the reported uncertainty is actually an underestimate of the correct one. 
We remind to the reader that this error come from the so called stability criterion, i. e. it is 
obtained looking at those approximants that minimize the difference between the estimates 
at the highest available orders (see Sec. II C). So, the considerable discrepancies between 
the four- and five-loop results lead to serious doubts in the existence of the cubic fixed point 
in two dimensions. 

In order to understand better these somewhat strange results, we report the values of 
the coordinate u* obtained for the cubic fixed point using several Pade approximants for 
iV = 3,4,8; the estimates are presented in Fig. 1 as functions of 6. Let us consider first 
the case N = 3 as a typical example. If one limits himself with only three lower-order 
approximants [2/1], [3/1], and [2/2], he easily finds that they minimize their differences 
under b ~ 0, leading in such a way to the estimate u* ~ 0.7. Taking into account two 
non-defective Pade approximants [4/1] and [3/2] (note the oscillating behavior of the [3/2] 
approximant for b < 1, signaling the presence of close singularities) existing at the five-loop 
level shifts the zone of stability to b ~ 2, thus leading to the estimate u* ~ 0.5. Moreover, 
the approximant [4/1] with 6 = 1, that is usually considered as one of the best approximants, 
results in the estimate u* = 0.050 very close to zero. Because of the alternative character of 
RG expansions, it looks very likely that the unknown six-loop contribution (and the higher- 
order ones) will locate the stability region somewhere near 6 ~ 1, leading finally to the 
coalescence of the Ising and cubic fixed point. According to this scenario, the cubic fixed 
point, found at finite order in perturbation theory, is probably only an artifact due to the 
finiteness of the perturbative series. 

The same scenario is possible also for other values of N. From Fig. 1 we see that the 
region of maximum stability always shifts from 6 ~ to 6 ~ 2 with increasing the order 
of approximation from four to five loops, moving the coordinate of the fixed point from 
u* ~ 0.5 for N = 4 (u* ~ 0.2 for JV = 8) to u* ~ 0.3 (u* ~ 0.1). Let us stress again that the 
value given by the approximant [4/1] with 6 = 1 is always very close to zero, as is seen from 
Table VI. Note that, with increasing N, the distance between the cubic fixed point and the 
Ising one reduces rapidly. 

In the limit A" — > oo the series simplify as at the four- loop level (see Ref. [3]). We only 
mention that with increasing the length of the RG series the coordinate u* of the cubic fixed 
point shifts from u* ~ 0.08 to u* ~ 0.03 that again is much more close to zero. 
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N=0 



N=3 




U *0.2 



N=4 
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0.6 
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-0.2 
-0.4 
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0.5 1 1.5 2 2.5 3 

b 

FIG. 1. Coordinate u* of the cubic fixed point within Pade-Borel method for some values of N. 

C. The random O(M) model 

The MN model in the limit iV — > describes the critical behavior of quenched dilute 
O(M) models. Being w proportional to minus the variance of the disorder [18,22] only the 
region with u < is of physical interest. We have already shown that for M > 2 the O(M) 
fixed point is stable. To analyze the RG flow in the whole physical region we do not use here 
advanced resummation procedures developed [30] to avoid Borel non-summability at fixed 
u/v [52], but limit ourselves by a simple Pade analysis, since it is sufficient for our aims. 
Within this method we check for several M > 2 that no other FP with u < exists when 
the resummation is effective. Thus, whenever the transition of dilute models is second-order 
(i.e. under the percolation threshold), the critical behavior is of O(M) type. An unstable 
FP is found in the unphysical region u, v > 0, which delimits the domain of attraction of the 
O(M) and O(0) fixed points, in agreement with the extrapolation of e-expansion at e = 2 
[17]. 

Finally, we discuss the fate of the random fixed point governing the critical behavior of 
the weakly-disordered Ising model (M = 1 and N = 0) that is found in three dimensions 
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[9,30,22] and close to four dimensions by means of the ^fi expansion [57,8,22]. In this case 
the Ising FP is marginally stable for u > [2], so the RG flow is driven to u* = from the 
physical region (in the unphysical region the Ising FP is marginally unstable and the RG 
flows run into the O(0) FP). We find, for the majority of the considered approximants, a fixed 
point with negative u*, as reported in Fig. 1. A possible estimate, according to stability 
criteria is u* = —0.1(1), but if we concentrate on some certain approximants we obtain 
u* = -0.090 for the [3/1] and u* = -0.030 for the [4/1] (both with 6 = 1). In particular, 
the last value is very close to zero, i. e. to the value predicted by the asymptotically exact 
solution that has been obtained in the framework of the fermionic representation [5,6,2]. It is 
worthy to note that the five-loop results for N = seem to be less scattered than analogous 
four- loop estimates obtained by means of Chisholm-Borel resummation technique [58], and 
they look more precise than their five-loop counterparts for finite N. 

D. The cubic model for N = 2 

The four-loop analysis of Ref. [3] for N = 2 turned out to be compatible with the presence 
of a line of the fixed points joining the 0(2)-symmetric and the decoupled Ising fixed points. 
The lines of zeros of the two j3 functions were found to be practically parallel and the quoted 
error was bigger than distance between them. This line of fixed points with continuously 
varying critical exponents is in agreement with what is expected from the correspondence, 
at the critical point, between the cubic model and the Ashkin- Teller and the planar model 
with fourth-order anisotropy [1,3]. We are now in a position to verify this statement at the 
five-loop level. 

First, we analyze the series with the conformal mapping method. Again we find that 
zeros of two (3 functions form two parallel lines, while the apparent uncertainties seem to be 
smaller than their separation. Of course, this fact may simply indicate that the model has 
no fixed point at all. Let us, however, look more carefully to this result and, in particular, 
to the accuracy of the quoted error. In fact, as we have already seen for N > 2, the error 
coming from stability criteria is likely an underestimate of the correct one. To understand 
better the situation, we use the Pade-Borel method. In Fig. 2 we report the curves of zeros 
of the two (5 functions given by several Pade approximants under 6=1, the value that for 
N > 2 always leads to good results and that is the best for the fixed point values of O(N) 
and Ising model [15]. The four approximants for fa are always well-defined. They are hardly 
distinguishable close to the 0(2) fixed point and their separation slowly increases moving 
toward the v axis. The coordinate of the Ising fixed point v ~ 1.8 is obtained using the [4/1] 
approximant, since approximants of [L-l/1] type proved to give rather precise estimates for 
the fixed point location both in two and three [59,60] dimensions. The situation is a bit 
worse for the function fa. In fact, the working approximants are well defined close to the 
Ising fixed point, but approaching the u axis they becomes defective. The approximant 
[3/2] starts oscillating around u ~ 0.8, while [4/1] is bad in the range u G [1, 1.5] and [3/1] 
for u > 1.3. Also the values of zeros of the approximant [4/1] for u > 1.5 are not reliable 
enough, since they may suffer of the effect of close singularities. 

Despite of these shortcomings, we can obtain a rich information from Fig. 2. Indeed, the 
line of zeros of fa given by the approximant [4/1] practically coincide with those of the fa 
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u 

FIG. 2. Zeros of (3u (continuous lines) and (3y (points) for several Pade approximant (all with 
6 = 1). 

from the Ising fixed point up to u ~ 0.8. For greater u various approximants for results 
in the lines of zeros that diverge leaving, however, the line [4/1] of j3u zeros between them. 
Keeping in mind a finite length of the RG series and the influence of non-analytic terms 
missed by the perturbation theory, we retain this fact as a strong evidence in favor of the 
continuous line of fixed points. The best estimate for this line is believed to be that given by 
the approximant [4/1] for Thus it will be used in what follows to calculate continuous 
varying critical exponents. 

We evaluate the smallest eigenvalues of the Q matrix along the line of fixed points both 
with conformal mapping and Pade-Borel method. We find that it is always compatible with 
zero (with the uncertainty of the resummation), that is a necessary condition for having a 
stable line of fixed points. 

When evaluating the critical exponents, one should keep in mind that the limit z — > is 
not simply accessible perturbatively since it corresponds to the two-dimensional XY model 
which is known to behave in a quite specific manner. In particular, this model does not 
undergo an ordinary transition into the ordered phase at any finite temperature and its 
critical behavior is essentially controlled by the vortex excitations [61]. Such excitations 
lead to an exponentially diverging correlation length at finite temperature that can not be 
accounted for within the A0 4 model Eq. (1.1) dealt with in this paper. Since arriving to 
the point z = a new physics emerges, it is natural to assume that the observables as 
functions of the fixed point location may be non-analytic near this point. Hence, what we 
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FIG. 3. The exponent y as function of x from free and constrained resummation. 
"first approximation" is Eq. (3.7) and the line "second approximation" is Eq. (3.8). 



The line 



can really explore trusting upon our five-loop expansions is a domain corresponding to finite 
(and not too small) values of z. Oppositely, the limit z — > oo (or 1/z — » 0) looks quite 
" undangerous" in the above sense since it corresponds to the critical behavior close to that 
of the Ising model, which was shown not to be influenced considerably by the non-analytic 
terms even in two dimensions [15]. Note that just presented results concerning the line of 
the fixed points confirm this idea. Indeed, as seen from Fig. 2, at the "Ising side", i. e. for 
< u ~ 0.8, the zeros of both (3 functions form smooth curves running very close to each 
other. The closer the " XY side", however, the stronger the estimates for (3 function zeros 
are scattered, indicating, likely, the increasing impact of non-analytic contributions. 

The expected value of r\ is 1/4 independently on the location of a fixed point within the 
line. The best way to check the constantness of r\ along the line of fixed points is probably 
to write the RG function in terms of u and s = u + v, i.e. r] s (s,u) = r)(u, s — u). Then 
one resum the difference A(s,u) = r] s (s,u) — i] s (s,0). Along the line for all the five-loop 
approximants we always find |A(s,w)| < 8 x 1CT 3 . This leads us to conclude that the two 
dimensional LGW approach is able to keep the constantness of r\ within an error of about 
3%. Note that the previous quoted problems concerning non-analyticities close to the 0(2) 
side do not significantly affect the estimates of r\. 

For the exponent y = rj — T) t = 1/(2 — v) it was conjectured in Ref. [3] that it should 
behaves like 

2 2 v* . . 

y = - — ; — , where, x = — arctan — . (3-7) 



1 + x' 



7T 
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A direct reliable quantitative estimate of this exponents is impossible because of the effect of 
non-analyticities, in particular close to the 0(2) fixed point for the reasons explained above. 
In fact, we know from Ref. [15,40] that the resummation of the series for y at the Ising fixed 
point provides y = 1.02, very close to the exact value 1. Instead, at the 0(2) fixed point one 
has y ~ 1.25, that is quite far from a diverging z/, i.e. y = 2. A direct resummation of this 
exponent is reported in Fig. 3, and as expected it seems to reproduces the correct critical 
behavior only close to the Ising fixed point. 

In Ref. [3] it was proposed to constrain the exponents to assume the exactly known value 
along the axes to have better quantitative results. We apply here a different method of 
constrained analysis with respect to Ref. [3] . We prefer to constrain the series of y expressed 
in terms of u and s = u+v (as previously for the exponent 77), since the results obtained with 
this constrain appear more stable. The constrained Pade-Borel approximants are reported in 
Fig. 3. At the Ising side, they are practically indistinguishable from the unconstrained ones 
up to x ~ 0.6, but then they start oscillating, signaling the presence of singularities leading to 
bad quantitative estimates. The conformal mapping results are practically equivalent. The 
numerical data thus obtained are not in agreement with Eq. (3.7), even if this reproduced 
the right qualitative behavior. A quadratic form in the denominator of y as e.g. 



fits the most of the data much better, as shown in Fig. 3. Unfortunately, we are not able 
to estimate the goodness of our resummation and, as a result, to verify Eq. (3.8). Perhaps, 
the exact behavior of the exponent v along the line of fixed points requires new method of 
analysis of the perturbative series and, in any case, it deserves for different studies on the 
subject like Monte Carlo simulation or high temperature expansion. 



The critical behavior of the MN model for M, N > 2 and v q > can be understood 
by means of non-perturbative arguments. In fact, as we explained in the Introduction, it 
turns out that the O(M) fixed point is always stable and oppositely the O(MN) one is 
always unstable. Thus the RG flow is driven to the O(M) fixed point from both positive 
and negative u. The only extra fixed point is provided by the e expansion with u* < 
which limits the attraction domain of the O(M) one. These features are reproduced by 
the resummation of the five-loop series. In particular, we found only one fixed point with 
non- vanishing coordinates that is the mixed one of the e expansion, that in fact turns out 
to be unstable. 

More interesting is the RG flow for v < 0. In this case the e expansion does not 
provide any fixed points, thus the transition (if any) cannot be continuous. Nevertheless 
for M = N = 2 we can relate the MN model to the 0(2) x 0(2) one, arguing that a fixed 
point with negative v* (located at u* ~ 4.6 and v* ~ —4.0 [36]) should exist, describing a 
finite temperature phase transition in the chiral universality class. Note that, as stressed 
in Refs. [34-36] (see also [51] for the three-dimensional analogous) such fixed point is in 
the region of non-Borel summability: in fact, z* = v* /u* ~ —0.86 < z± — —0.833, where 



4 



(3.8) 




E. The MN model for N, M > 2 
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z\ = Rmn/Rn delimits the region where the series are Borel summable (see Sec. II C). 
Anyway, in the course of the resummation with the conformal mapping, the singularity of 
the Borel transform closest to the origin has been taken into account. The resulting sequence 
of approximations is thus only asymptotic (as for the random models), but it is expected 
to provide a reasonable estimate as long as the resummation point is far from the region 
Z2 = z±2N/(N + 1), where the singularity on the real positive axis becomes the closest to 
the origin. 

With this in mind, we resum the (3 functions for several values of M and N and search for 
new fixed points with v* < 0. We consider for each (3 functions the 18 conformal mapping 
approximants with a — 0, 1, 2 and b = 5,7, 9, 11, 13, 15 that (as for the O(M) x O(N) model 
[36]) turned out to be the more stable with changing the number of terms considered in 
perturbation theory. Up to three loops we do not find a fixed point for any value of N and M, 
but at four and five loops a fixed point with z* very close to z\ appears for the majority of the 
324 couples of approximants of the (3 functions. In Table VII we report the five-loop results 
for the location of the fixed point u*,v* and the percentage of the 324 approximant having 
it. In the Table are reported also the values of z* (obtained independently in the course of 
the averaging procedure) and the values of z\ and z<i to make the comparison between them 
clear at first sight. The reported estimates are the averages over those approximants having 
a fixed point and the error bars are the variances. 

For M = N = 2 we reproduce the result of the chiral model [36]. Several features for 
other value of M and N can be extracted from Table VII. For each M there is a maximum 
N, we call N*(M), for which the FP exists. This is true up to a maximum value of M, 
we call M* , after that the results of the e expansion and large M are recovered. A reliable 
determination of M* and N*(M) is difficult. In fact, one has arbitrarily to fix a confidence 
level of the percentage for which the existence of the fixed point can be considered firm. 
For example, we can decide that the fixed point is credible when the percentage of the 
approximants displaying it is greater than the 70% and its existence is improbable when 
it is less than 40%. Within this method we have 2 < JV*(2) < 4, 3 < JV*(3) < 5, etc. 
Regarding M* the situation is even more difficult. With increasing M, the FP moves closer 
and closer to z<i, and, in fact, its existence is not really reliable. Anyway, we retain quite 
safe to state M* < 22, but probably its value is much lower. 

We check the stability of the results at four and five loops. For small values of M the 
location of the FP is rather stable, but it gets worse with increasing M, signaling that the 
estimate in this case is not so robust. 

We evaluated the eigenvalues of the Q matrix for each couple of approximants at their 
common zero. They always have a quite large positive real part, so the FP is stable. However, 
it is very hard to obtain reliable numerical estimates since the values strongly oscillate. At 
the same time, the focus behavior seems to be a peculiarity of the chiral model M = 
N = 2; for larger M, N only few approximants lead to a non-vanishing imaginary part. 
Even the estimates of the critical exponents are practically impossible probably because of 
non-analyticities and the large values of the coupling constants involved in the calculation. 
Anyway, the exponent 77 for moderate M seems to be compatible with 1/4 that is the value 
of the Ising and XY models. 

Let us comment that a parallel analysis in d = 3 [28] seems to indicate a similar structure 
of fixed points, but the results are more stable (as usual) both because of the knowledge of 
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TABLE VII. Location (u*,v*) of the stable fixed point with v* < at five loops for several M 
and N. The estimate of z* = v* /u* and the value of z\ and Z2 (see the text) are also reported. 



M 


N 


% 


u* 


V* 


z* 




Z2 


2 


2 


82 


4.65(18) 


-3.93(33) 


-0.85(5) 


-0.83 


-1.11 


2 


3 


57 


4.66(41) 


-3.73(46) 


-0.80(5) 


-0.71 


-1.07 


2 


4 


42 


4.88(49) 


-3.72(51) 


-0.76(5) 


-0.62 


-1 


2 


5 


31 


4.99(51) 


-3.64(44) 


-0.73(5) 


-0.56 


-0.93 


3 


2 


86 


5.37(25) 


-4.34(29) 


-0.81(4) 


-0.79 


-1.05 


3 


3 


65 


5.20(60) 


-3.89(53) 


-0.75(5) 


-0.65 


-0.97 


3 


4 


44 


5.34(62) 


-3.79(54) 


-0.71(5) 


-0.55 


-0.88 


3 


5 


36 


5.81(82) 


-3.95(65) 


-0.68(5) 


-0.48 


-0.80 


4 


2 


89 


6.04(39) 


-4.71(26) 


-0.78(4) 


-0.75 


-1 


4 


3 


67 


5.64(79) 


-4.01(61) 


-0.71(5) 


-0.6 


-0.9 


4 


4 


52 


6.03(97) 


-4.05(76) 


-0.67(5) 


-0.5 


-0.8 


5 


2 


91 


6.63(57) 


-5.05(31) 


-0.76(4) 


-0.72 


-0.96 


5 


3 


67 


5.94(89) 


-4.07(58) 


-0.69(5) 


-0.56 


-0.85 


5 


4 


50 


6.5(1.2) 


-4.15(77) 


-0.64(4) 


-0.46 


-0.74 


6 


2 


90 


7.10(75) 


-5.32(40) 


-0.75(4) 


-0.7 


-0.93 


6 


3 


67 


6.2(1.0) 


-4.13(61) 


-0.67(5) 


-0.53 


-0.81 


7 


2 


93 


7.5(1.0) 


-5.55(51) 


-0.74(4) 


-0.68 


-0.91 


8 


2 


88 


7.7(1.1) 


-5.68(66) 


-0.74(3) 


-0.67 


-0.89 


8 


3 


64 


6.5(1.2) 


-4.2(6) 


-0.58(2) 


-0.5 


-0.75 


9 


2 


84 


7.9(1.2) 


-5.83(70) 


-0.74(3) 


-0.65 


-0.87 


10 


2 


83 


8.0(1.2) 


-5.87(70) 


-0.74(3) 


-0.64 


-0.86 


15 


2 


75 


8.1(1.0) 


-6.01(55) 


-0.74(3) 


-0.61 


-0.81 


15 


4 


15 


10(2) 


-5.2(9) 


-0.52(1) 


-0.34 


-0.54 


20 


2 


55 


8.9(1.0) 


-6.57(54) 


-0.74(2) 


-0.58 


-0.78 


20 


3 


55 


8(1) 


-4.6(5) 


-0.58(2) 


-0.41 


-0.62 


22 


2 


23 


9.2(5) 


-6.7(3) 


-0.73(2) 


-0.58 


-0.77 


25 


2 


21 


10(1) 


-7.5(5) 


-0.73(2) 


-0.57 


-0.76 



six-loop terms [30] and because of the weaker effect of non-analyticities (in particular, for 
the estimates of the exponents). 

To conclude this section, we want to stress that Table VII should be read with care. The 
point is that all the found FP have z* < z\ being located in the region where the series are 
not Borel summable and all the FP with z* close to Z2 are not credible because of the bad 
behavior of the resummed approximants close to this region. 

IV. CONCLUSION 

To summarize, the critical behaviors of the two-dimensional iV-vector cubic and MN 
models have been studied within the renormalizat ion-group approach. The five-loop con- 
tributions to the P functions and critical exponents have been calculated and the five-loop 
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RG series have been resummed by means of Pade-Borel-Leroy procedure and the conformal 
mapping technique. 

For the cubic planar model (N = 2) we have found that the continuous line of fixed points 
connecting the Heisenberg and the Ising ones is well reproduced by the resummed five-loop 
RG series. Moreover, the five-loop terms being taken into account make the lines of zeros of 
f3 functions for u and v closer to each another thus improving the results of the lower-order 
approximation. For the cubic model with N > 2, the five-loop contributions have been 
shown to shift the cubic fixed point, given by the four-loop approximation, towards the 
Ising one. This may be considered as an argument in favor of the idea that the existence of 
cubic fixed point in two dimensions for N > 3 is an artifact of the perturbative analysis. 

The models with N = describing the critical thermodynamics of two-dimensional 
weakly-disordered O(M) systems has been also studied. The results obtained have been 
found to be compatible with the conclusion that in two dimensions the impure critical 
behavior is governed by the O(M) fixed point, even in the Ising case, where ati = and the 
Harris criterion is inconclusive. 

The five-loop RG analysis of the two-dimensional MA-model with M, N > 2 has been 
also performed. For v > we reproduced all the non-perturbative results. It was shown, 
in particular, that the transition is driven to the 0(M) fixed point, that only for M < 2 
describes a finite temperature phase transition. We also found a stable fixed point in the 
region with v q < that has no counterpart in e and large M expansion, but its location 
is in the region of non-Borel summability of the series and its existence is still doubtful. 
At fixed M, this new fixed point is found for N < N*(M) up to a maximum value M*, 
after which it disappears. Whether this fixed point describes a finite temperature phase 
transition (allowed because of the discrete symmetry Cn) or a zero temperature one cannot 
be discerned by our analysis; only lattice techniques such as Monte Carlo simulation and 
high temperature expansion or real experiments can completely clarify this point. 
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